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In a recent paper [Phys. Rev. A 91 , 053844 (2015)], Mukamel and Dorfman compare spectroscopies 
performed with classical vs. quantum light, and conclude that nonlinear quantum-spectroscopy sig¬ 
nals cannot be obtained from averaging their classical-spectroscopy counterparts over the Glauber- 
Sudarshan quasiprobability distribution of the quantum field. In this Comment, we show that this 
interpretation is correct only if one assumes that a classical signal is given by a classical approxi¬ 
mation for the held. While such an assumption can be useful for comparing theoretical results, it 
is never realized in laser spectroscopy experiments that typically use coherent states. Thus, instead 
of using classical signals, the connection between coherent states and quantum states of light must 
be considered. We rigorously show that quantum spectroscopy can always be projected from the 
experimentally realized coherent-state spectroscopy regardless how nonlinear the system response 
is. 


Finding and understanding the border between the 
classical and quantum regimes is essential for establishing 
the foundations of modern optics. For example, the sem¬ 
inal contributions by Glauber[lj and SudarshanQ iden¬ 
tified coherent states, |/3), as the most classical form of 
quantized light. Once the phase of a single-mode light 
field, such as is emitted by a perfect laser, is defined 
against a reference (otherwise light’s phase remains ran¬ 
dom, as pointed out by MplmerQ), it is quantum sta¬ 
tistically described by one coherent state. Normally- 
ordered expectation values for \(3) are obtained by a clas¬ 
sical factorization [3, where the photon operators B 
and are replaced by complex amplitudes /? and /?*, 
respectively, as if the field had a definite classical phase 
and amplitude. However, coherent states have quantum 
fluctuations that invalidate the classical factorization for 
expectation values that are not normally ordered. 

One can use coherent states to express the density ma¬ 
trix of a single-mode light via[ll Q 

p = j cPpP{P)\p){p\, (1) 

where F(/?) is the Glauber-Sudarshan function. This re¬ 
lation is the starting point of the analysis by Mukamel 
and DorfmanQ. As pointed out by Glauber and Su- 
darshan, Eq. © appears to be a probabilistic phase- 
space average over coherent states, i.e. |/3)(/3|. For a 
light described by |/3o), P{li) becomes 5{j3 — /3o). States 
of light with a positive-definite P(/3) must necessarily 
have equal or greater fluctuations compared to coherent 
states. Since their fluctuations result from a probabilis¬ 
tic average of the coherent states, these states of light 
are often referred to as “semiclassical“. There are also 
many states of light, such as Fock and Schrodinger-cat 
states, whose F(/?) is partially negative valued or even 
nonanalytic. Such sources are often referred to as “true 


quantum light". 

At the moment, there exist very few experimental re¬ 
alizations of true quantum light, which are technically 
challenging and usually limited to very low intensities. 
Therefore, it is interesting to determine whether spec¬ 
troscopy performed with semiclassical vs. true quantum 
sources are connected. In this context, experiments per¬ 
formed with coherent states have been labeled as “clas¬ 
sical spectroscopy”, whereas those performed with true 
quantum light have been labeled “quantum (optical) 
spectroscopy” [3 0|. However, it is important to bear 
in mind that since “classical spectroscopy” refers to mea¬ 
surements made with coherent states, there are still quan¬ 
tum aspects to the measurement, namely the intrinsic 
quantum fluctuations, and the results of these measure¬ 
ments are not necessarily the same as for using the ap¬ 
proximation of purely classical light, which is actually not 
experimentally possible. We will try to emphasize this 
point by avoiding the term “classical spectroscopy” and 
explicitly stating when we are referring measurements 
made with coherent states. 

In his 1963 paper, Sudarshan essentially writes a pro¬ 
jection between outcomes for coherent-states and quan¬ 
tum light [ 2 ]. This projection can be used to establish a 
connection between quantum spectroscopy and coherent- 
state spectroscopy. Using the notation from Refs. [i-ll|, 
the system response to quantum light, Rqm, in terms of 
the responses to coherent-states, R\p), is 


Rqm = J d^pPWRip), 


( 2 ) 


where the quantum-light statistics are defined by P{(3). 
Sudershan defines via normally ordered operators, 
and we show later that, in general, the outcome of a spec¬ 
troscopy measurement can eventually be expressed via 
R\/ 3 )- In the case where a spectroscopic technique uses 
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multiple single-mode sources or a multi-mode source, the 
connection between coherent-state and quantum-optical 
measurements, m, can be generalized by replacing /3 by 
/3 i,/ 32 etc., where f3j defines the quantum statistics of 
each orthogonal mode applied in the technique. We will 
next study whether the claims made by Mukamel and 
Dorfmanp that Eq. (I2D “only holds for the linear re¬ 
sponse, and does not apply to the nonlinear response” 
and “The quantum response carries additional valuable 
information about response and spontaneous fluctuations 
of matter that may not be retrieved from the classical re¬ 
sponse by simple data processing” are justified. 

In their paper, Mukamel-Dorfmanj^ consider two ex¬ 
amples where i?|^) is first computed and then it is ana¬ 
lyzed how the projection ([2]) connects different levels of 
calculations. In the first example, they discuss the third- 
order nonlinear spectrum of a system which they denote 
as S{uj). They apply a classical factorization, i.e. they 
replace the photon operators by complex-valued ampli¬ 
tudes, to obtain which they treat as the system 

response to a classical field in Eq. (14) of Ref. Q. They 
compare this result with a calculation that does not apply 
the classical factorization, giving the quantum response 
S{uj)qm- They correctly find that Sqm cannot be repro¬ 
duced by the projection ([2]) when i?|^) is replaced by Sd- 
Furthermore, they argue that this violation results be¬ 
cause “a classical field does not change in the course of 
field-matter interactions”. 

To analyze this apparent violation of the quantum- 
spectroscopy concept 0, it is helpful to consider 

whether a laser-spectroscopy experiment performed with 
a coherent state, actually yields Sd {Rd in the notation 
we are using) resulting from a classical factorization. In 
this context, it is a very instructive example is to consider 
the laser excitation of a two-level system inside a cavity. 
Using the Jaynes-Cummings model[l2| to compute the 
system response to a coherent state |/3) without the 
classical factorization shows that the occupation of the 
excited state first increases but then collapses [l^ even 
without damping. It was also found in 1960s that the 
occupation collapse is followed by a series of revivals and 
further collapses [l^. However, as soon as a classical 
factorization is applied, the theory reduces to the optical 


Bloch equations!^ Il5j that produce the occupation Rd 
corresponding to undamped Rabi oscillations when de¬ 
phasing is excluded. Consequently, the response, to 
a coherent state is drastically different from Rd obtained 
with a classical factorization. Obviously, the usefulness of 
vs. Rd can be judged based on which of them pre¬ 
dicts the outcome of an actual laser-spectroscopy mea¬ 
surement. In 1987, Rempe et al. 0 used a maser to 
measure a clear sequence of collapses and revivals, which 
demonstrates that a maser/laser spectroscopy measures 
R\I3), not Rd- 

This conclusion is not surprising because a measure¬ 
ment cannot switch off the quantum aspects of a coher¬ 


ent state — at its best, a laser is in maximally classical 
coherent state |/3) whose quantum fluctuations still influ¬ 
ence the measured response R|/ 3 ). As originally formu¬ 
lated by SudarshanQ, the projection relation ([2]) con¬ 
nects R|/ 3 ), not Rd, with a quantum-optical response, 
Rqm, to a true quantum-light source. In our quantum- 
spectroscopy investigations [sMllj. we carefully apply this 
fundamental connection, and correctly use i?|^) as the 
system response to a coherent state, obtained without 
any classical approximations to the field and response. 

By using classically approximated Rd instead of the 
full R\p), Mukamel and DorfmanQ conclude that the 
projection © “breaks down in the nonlinear regime”. 
Therefore, it is instructive to check whether this conclu¬ 
sion holds for spectroscopy using coherent states. Before 
the light excites matter, it can be assumed not to be 
entangled with the matter such that the initial state of 
the total system can be expressed as ptot (0) = Pmat ® p 
where Pmat is the initial state of matter while p defines 
the light state. For laser spectroscopy, p becomes |/3)(/3| 
of a coherent state, which yields /3|p)(0) = pmat ® |/3)(/3| 
as the initial state of the combined matter-light system. 
When this result and Eq. o are merged, the initial state 
becomes 

/5tot(0) = y'd2/3P(/3)p|p)(0). (3) 

For a closed system consisting of only the matter and 
the exciting light, the full evolution follows from the 
system Hamiltonian H including the quantized light- 
matter coupling as well as possible many-body inter¬ 
actions within the matter. In the Schrodinger picture, 
the exact time evolution formally follows from ptotit) = 
/5tot(0) regardless how nonlinear the in¬ 

teractions are. By inserting Eq. (|31) into this relation, we 
find 

Ptot{t) = j (fpP{P) p\p){t) (4) 


because P{a) is not an operator. In this decomposition, 
P\p){t) corresponds to the exact evolution following an 
excitation by a coherent-state \f3). It is clear that the 
part of P|/ 3 )(t) describing the light does not necessarily 
remains a pure \I3) when interactions are present, which 
invalidates the classical factorization. Therefore, o is 
no longer valid when the classical factorization is applied 
to approximate the exact p\p){t)- 

The general system response Rqm corresponds to an 
expectation value of an observable R, i.e. Rqm = (R) = 


Tr 


R Ptot (t) 


, or a Fourier transform over t yielding spec¬ 


tral information. By using the connection (j4]), we find 


Rqm = Tr 


R / d2^R(^)p|p)(t) 


= / cPpP{f5)TY Rpip)(t) 


(5) 
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which follows by exchanging the trace an integration. 


is the exact system 
’). Therefore, Eq. JS]) 


We notice now that Tr i?p|,g)(t) 
response to a coherent state |/ 
reduces to an exact projection ([2]) even when the spectro¬ 
scopic response is nonlinear. The opposite conclusion by 
Mukamel and Dorfmanj^ stems from the classical factor¬ 
ization which, however, is not applicable to describe real 
laser spectroscopy. We have shown in Ref. Q that the 
projection ([2]) is also valid for open systems. Basically, 
the projection is a direct consequence of the linearity of 
quantum mechanics, irregardless of the nonlinearity of 
R\p) induced by many-body and quantum-optical inter¬ 
actions. 

In their second example, Mukamel and DorfmanQ in¬ 
vestigate nonlinear fluctuation-dissipation relations with 
superoperators. In their Eq. (19), the system response 
is Tr [Oip[t)\ where Oi correspond to our R and p{t) to 
Ptot{i) converted to the interaction picture. They ap¬ 
ply a classical factorization when computing Tr [Oip{t)\ 
and conclude that “The deep reason why the CRF and 
QRF are not simply related is in the lack of a fluctuation- 
dissipation relation in the nonlinear regime” where CRF 
(QRF) stands for the classical (quantum) response func¬ 
tion. However, a laser spectroscopy experiment does not 
measure Tr [Oip(t)] obtained with a classical factoriza¬ 


tion for the field, but rather the full Tr 


Rpmit) 


re¬ 


sponse to a coherent state. As we show above, the pro¬ 
jection ([2]) connects laser spectroscopy performed with 
coherent states and quantum spectroscopy. 

Since experiments realizing laser spectroscopy with co¬ 
herent states are common place whereas the technol¬ 
ogy for generating arbitrary states of quantum light is 
still primitive, the projection m yields an interesting 
method for deducing quantum responses from a large 
set of laser-spectroscopy measurements using coherent 
states. However, many interesting quantum sources have 
an extremely complicated /*(/?) function consisting, e.g., 
of high-order derivatives of a delta function. Since 
cannot be measured experimentally throughout the en¬ 
tire phase-space and any measurement inevitably con¬ 
tains some additional noise, it is clear that integrating 
@ with experimental data and a rapidly changing (or 
nonanalytic) quantum P(/3) is challenging. 

We have devised a numerical algorithm [3, based on a 
cluster-expansion transformation (CET)[l7|, to robustly 
integrate f (Pp P{P) R\p) from a finite number of |/3) 
measurements. We have thoroughly tested our GET 
projection in Ref. Q with the Jaynes-Cummings model, 
which is one of the few nontrivial examples where R\p) 
can be computed exactly. More specifically, we computed 
R\p) for 240 P values for 0 < |/3| < 9 resulting to the 
well-known collapse-revival sequences explained above. 
We then applied our GET projection to determine Rqm 
to Fock-states as quantum sources and fully recovered 
the related quantum-Rabi oscillations, which also have 


been measured[l3. In the same context, we have also 
shownQ that the GET projection method works even 
when random noise is added to the data demonstrating 
that our method is robust against experimental noise. 
Hence, it can be directly applied to project coherent- 
state laser spectroscopy measurements to true quantum- 
optical spectroscopy. 

We have also shown in Ref. Q that the quantum 
and classical spectroscopies produce different results only 
when exhibits nonlinearities. Therefore, the dif¬ 
ference of i?QM and R\p'^ is a particularly sensitive 
method for detecting many-body and quantum-optical 
effects that produce nonlinearities. We have applied 
the GET projection to detect complex many-body quasi- 
particles such as biexcitonsQ and dropletons[l3 in a 
semiconductor quantum well as well as a photon-density 
correlation [llj in the semiconductor quantum-dot emis¬ 
sion, and we expect our method to be an extremely use¬ 
ful extension of laser spectroscopy into quantum spec¬ 
troscopy. 

In conclusion, Mukamel and DorfmanQ analyze how 
different leyels of theory compare with one another, and 
we fully agree that a computation performed with a 
classical factorization does not contain the quantum- 
statistical information needed to project quantum-optical 
responses. Howeyer, a laser-spectroscopic measurement 
cannot be described by the classical factorization since it 
cannot switch off the quantum fluctuations of light, such 
that it always contains the quantum-statistical informa¬ 
tion required to yalidate the projection ([2]). We haye 
shown that laser spectroscopy performed using coher¬ 
ent states connects with measurements performed with 
true quantum sources yia the projection ([2|), which is ex¬ 
tremely useful because one can routinely realize coherent- 
state spectroscopy while quantum- source deyelopment is 
challenging. Eyen though quantum spectra can be pro¬ 
jected from a massiye set of coherent-state spectroscopy 
data, it is not any less exciting. One may view that 
a set of laser excitations disperse the full many-body 
and quantum-optical information between all possibili¬ 
ties, whereas quantum spectroscopy condenses this in¬ 
formation to an isolated quantum response. Analogous 
“parallelization” of outcomes makes quantum computing 
interesting compared to its classical counterparts. 

We thank Shaul Mukamel for discussions and providing 
us a preprint of Ref. The work in Marburg is sup¬ 
ported by the Deutsche Forschungsgemeinschaft (grant 
KI 917/2-2). 
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